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ABSTRACT 

Lensing flux-ratio anomalies have been frequently observed and taken as evidence for 
the presence of abundant dark matter substructures in lensing galaxies, as predicted 
by the cold dark matter (CDM) model of cosmogony. In previous work, we exam- 
ined the cusp-caustic relations of the multiple images of background quasars lensed by 
galaxy-scale dark matter haloes, using a suite of high-resolution A'^-body simulations 
(the Aquarius simulations). In this work, we extend our previous calculations to incor- 
porate both the baryonic and diffuse dark components in lensing haloes. We include in 
each lensing simulation: (1) a satellite galaxy population derived from a semi-analytic 
model applied to the Aquarius haloes, (2) an empirical Milky- Way globular cluster 
population and (3) satellite streams (diffuse dark component) identified in the sim- 
ulations. Accounting for these extra components, we confirm our earlier conclusion 
that the abundance of intrinsic substructures (dark or bright, bound or diffuse) is not 
sufficient to explain the observed frequency of cusp-caustic violations in the CLASS 
survey. We conclude that the observed effect could be the result of the small number 
statistics of CLASS, or intergalactic haloes along the line of sight acting as additional 
sources of lensing flux anomalies. Another possibility is that this discrepancy signals 
a failure of the CDM model. 

Key words: Gravitational lensing - dark matter - galaxies: ellipticals - galaxies: 
formation 



1 INTRODUCTION 

The cold dark matter (CDM) cosmogony predicts that cos- 
mic structures form hierarchically through a succession of 
mergers and accretions. The mass of the Milky Way's own 
dark halo is curr ently constrained to ~ 1 — 5 x 10^^ Mq (e.g. 
IGuo et al.|[2009l ). Numerical simulations predict a wealth of 
dark matter substructures (self-bound subhaloes) suryivin; 
in haloes of this mass (e.g.lGao et al.ll2004l : [Cao et al.|[20o". 
iDiemand et al.ll2008l : ISpringel et al.ll2008l '). These substruo 
tures have a power-law mass function; scaling this to a satel- 
lite galaxy luminosity function (LF) by adopting a fixed 
mass-to-light ratio consistent with the brightest objects 
overpredicts the number of satellite galax;ies of the Milky 
Way (MW) and M31 by a factor of sev eral hundred — the 
so-called 'missing satellite' problem (e.g. iKlvpin et al.lll999l : 
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I Moore et al. 1 1 19991 : for a recent review see lKravtsov|[2oTol '). 

For many years, theoretical models of galaxy formation have 
predicted that the star formation efhciency in low-mass 
haloes can be str ongly suppressed through a combination of 
photoionization ('E fstathioulll992l ) and supernova feedback 
(Whi te fc Rees 197^. This suppression renders many such 
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cent discoveries of low luminosity satellites i n the MW and 
M31 and corrections for in completeness (e.g. IKoposov et al.l 
l2008l : Irbllerud et al.ll200^ 1 imply reasonable agreement be- 
tween observations and theoretical predictions for MW-like 

I ' 11 1 

haloes IjBenson et a l. 2002). However, because they rely on 

uncertain baryonic physics and currently limited data, these 
comparisons so far only provide an indirect and uncertain 
test of the abundance and mass function of CDM substruc- 
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tures. Since gravitational lensing traces mass, regardless of 
its luminosity, it can provide a powerful alternative tool to 
investigate otherwise ' hidden' subhalqes and directly test 
the C DM model (e.g. iMacLeod et al.l l2009l : IVegetti et al] 
I2OO9I ). 

Flux-ratio anomalies observed in multiple images 
of lensed quasars are often cited as evidence for sub- 
struc tu res in lensing g a laxies (e.g. Mao fc Schneider] 
199i; iMetcalf fc Madaul I2001I: iMetcalf fc Zhaol 1200^^ 



Dalai fc Kochanekl pOO^: IChibal I2OO2I: iKochanek fc Dala 
2004; Metcalf et al.' '2004*: 'Chiba et al.' '2005!: ISugai et al 
2007; McKcan ct al. 2007; More ct al. 2009). S ome s tudies, 
e.g. jPalal fc Kochanekl (20021 ). iBradac et al.1 (j2004h using 
CDM simulations have concluded that the substructure 
population expected for a galaxy-scale lens can reproduce 
the observed fl u x anomaly cases. However , other studies, 
e.g. iMao et al.1 (120041') . lArnara et al.1 (|2006l '). [Maccio ct al,^ 
l|2006t) iMaccio fc Mirandal (l2006l ) "and our own previous 
work l|Xu et al.ll2009l l have argued that the simulated CDM 
substructure population is not sufficient to explain the 
observed frequency of lensing flux-ratio anomalies. This 
negative conclusion does not yet rule out the predictions of 
CDM, as the observed sample of 'anomalous' lenses is small. 
Intergalactic haloes along the line of sight may also perturb 
the l ensing poten t ial an d cause the observed lensin g anoma- 
lies (IChen et al.1 l2003l: IWambsganss et al.l l2005l: iMetcaU 



2005allbl: iMiranda fc Maccid l2007l: iPuchwein fc Hilbert 



20091 ). 



One extreme case 
cusp-caustic violation, 
five lens sys t ems of 
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of flux-ratio anomalies is the 
which has been observed in 
a cusp geometry: B0712-I-472 

B1422+231 (jlmpev et al.1 

Patnaik fc Narasimhal l200ll }. B2045-I-265 
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and RXJ0911-I-0551 (Bade et al. 1997; Burud ct alT ll998 
In the first three cases (those detected in radio observations) 
these have been attributed to substructure lensing, while 
the remaining t wo ar e thought to b e due to mi crolensing 
(Anguita ct alj l2008l : iMorgan et al] bOOd "). In IXu et al.1 
1 2009 ) (which we summarise below) we found that even 
with the well-resolved subhalo population from the Aquar- 
ius simulations, the observed cusp lenses still violate the 
cusp-caustic relation more frequently than expected for 
a concordance CDM model. However, our earlier work 
did not consider the effect of baryons concentrated at the 
centre of subhaloes, which may modify their density profiles 
and so alter their lensing properties. In this work, we 
carry out a more thorough analysis of the lensing effect of 
baryons in substructures, proceeding as follows. In Section 
2 we use a semi-analytic model of galaxy formation to 
predict the baryon content of subhaloes in each Aquarius 
simulation: we model the density distribution of baryons in 
each subhalo and re-run our lensing simulations to study 
cusp-caustic violations. In Section 3, to estimate the effects 
of possible baryonic concentrations below the resolution 
limit of the simulation, we further include an empirical 
Galactic globular cluster population (Harris 1996) in our 
model. Finally, as typical surface density fluctuations of 
only a few per cent would be suffici e nt to cau s e flux -ratio 
anomalies (|Mao fc Schneideil Il998l : iLi et al.1 I2OO6I ). we 
consider irregularities in the surface density of the main 
halo due to the debris of tidally disrupted satellites (the 



diffuse component) in Section 4. Section 5 provides a short 
summary and a brief discussion of our results. 



2 SEMI-ANALYTIC GALAXIES AND THEIR 
CUSP-CAUSTIC VIOLATIONS 

2.1 The Aquarius Simulations 

The Aquarius project (|Springel et al.l [2OO8I : ISpringel et al.l 
I2OO8I ) is a suite of coUisionless ACDM A'^-body simulations 
run with the GADGET-3 code. Each of the six simulations in 
the suite (labelled A-F) centres on a single halo of mass ^ 
lO^^/i"^M0. These haloes were selected for high-resolution 
multi-mass resimulation from a lOO/i"^ Mpc^ simulation 
with a uniform mass resolution: a subsequent high-resolution 
study of this parent volume (the Millennium II simulation) 
confirms that in almost all respects the Aquarius sa mple is 
representative of haloes of t his m ass (Bovlan-Kolchi n et al.l 
I2OO9I : iBovlan-Kolchin et al.ll2009l ') . 

All six haloes were simulated at a resolution correspond- 
ing to a particle mass of ~ 10*/i~^Mq and a softening length 
of ~ 50/i~^ pc, resulting in ~ 130 million particles within the 
virial radius of the main halo at z = 0. This resolution of the 
Aquarius project is referred to as 'level 2' and the six haloes 
at this resolution are labelled Aq-A-2, Aq-B-2, Aq-C-2, Aq- 
D-2, Aq-E-2 and Aq-F-2. At this resolution, subhaloes are 
resolved to a mass limit of mgub > IQ^ Mq (correspond- 
ing to bound associations of 20 particles identified by the 
SUBFIND algorithm). 'We use these level 2 simulations in this 
work. Simulations of one halo (Aq-A) at higher and lower 
resolutions have demonstrated that the abundance, radial 
distribution and internal structure of subhaloes relevant to 
our lensing calculations are numerically converged to high 
accuracy at the resolution of the level 2 simulations (full 
details of thes e exte nsive convergence studies are given by 
ISpringel et al] (|2008D and lNavarro et al.1 l|2010l )'). 

Throughout this work, the cosmology of our lensing sim- 
ulations is the same as that used for the Aquarius project, 
with a matter density f2m = 0.25, cosmological constant 
Q.A = 0.75, Hubble constant h = //o/(100kms"^ Mpc"^) = 
0.73 and linear fluctuation amplitude ag = 0.9. 



2.2 Halo lensing model 

In lXu et al.l (|2009| ). we studied cusp-caustic violation prob- 
abilities using the six galaxy-scale CDM haloes provided by 
the Aquarius simulations. Mock lensing signals from back- 
ground quasars were used to determine the observed fre- 
quency of fiux-ratio anomalies due to dark matter substruc- 
tures in these haloes. The mass resolution of Aquarius is 
at least two orders of magnitude higher than that of pre- 
vious simulations used for lensing studies, resolving bound 
subhaloes down to a mass of lQ^h~^ Mq. In addition to 
the simulated dark matter distribution, we added a Hern- 
quist profile 'galaxy' at the centre of each main halo. 'We 
renormalized the dark matter particle mass and adiabati- 
cally contracted the particle distribution t o account for this 



extra contribution to the potential (e . g. Barnes fc White 
1984 : iBlumenthal et aL I ll986l : iMo et al.1 1 19981 : ICnedin et al 



20041 ) 



A particle-mesh (PM) code was applied to calculate 
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Figure 1. The masses and radial distributions of semi-analytic galaxies (bl ue) and dark matter subhaloes (black) from Aq-A-2, as well 
as the observed Milky- Way globular clusters (red; using the catalogue from lHarri£lll996l) . The panel on the left shows the stellar mass 
Mbaryon S-S a function of the total mass Mtot^li which in the case of semi-analytic galaxies is the total mass of a satellite galaxy and its 
host subhalo. The dynamical masses of globular clusters come from their observed V-band luminosities, assuming a mass-to-light ratio 
to be 3. Such masses are equal to their stellar masses Aftaryom assuming a zero dark matter content. In general, globular clusters have 
lower dynamical masses than satellite galaxies, but the latter have higher mass-to-light ratios than the former. The middle and right 
panels present the projected and 3D radial distributions of the three components. Globular clusters are much more concentrated towards 
the galactic centre than satellite galaxies. At the Einstein radius of the main halo (2 ~ 4h~^ kpc), the number of projected globular 
clusters is about 1 ~ 2 times that of dark matter subhaloes. Satellite galaxies are rare in the projected central region. 



lensing potentials, deflection angles and magnifications from 
the simulated lensing hal oes. The numeric al accuracy of this 
code is described fully in IXu et al.1 (|2009| ). Near the critical 
curves, high magnification makes changes in the image prop- 
erties sensitive to small perturbations in surface density. To 
estimate the effects of numerical noise in these regions, we 
used our code to calculate the lensing properties of isother- 
mal ellipsoids generated by Monte Carlo sampling, using an 
equivalent number of particles to the level 2 Aquarius simu- 
lations (~ 10* particles within the radius r2ocEl). The numer- 
ical accuracy of the defiection angle, convergence (surface 
density) and magnification fo r these Monte C arlo realisa- 
tions are shown in Figure 4 of IXu et al.1 (I2009D . The uncer- 
tainties around the Einstein radius (at about 0.02 r2oo) are 
0.2% for the deflection angle, 1% for the convergen c e, and 
< 10% for the magnification. Figure 6 of IXu et al.l (|2009l ) 
presents the probability distribution function of the image 
flux ratio -Rcusp (Eq.[T]) for cusp sources with an image open- 
ing angle smaller than 90° in the 10*-particle Monte Carlo 
realisations. This distribution was found to be broader than 
the equivalent analytical calculation as the result of discrete- 
ness noise, inherent in mapping an A^-body density field to 
a mesh. 

To correct for this sampling noise, each main 
Aquarius halo (including a central 'galaxy') was mod- 
elled as an isothermal ellipsoid (e.g. iRusin et al] 



eiied as an isotnerma i einp s oid (e.g. iKusm et ai.l 
20031: iRusin fc Kochane^ l2005l : iKoopmans et al.1 l200a 



Gavazzi et alj l2007t) . the lensing properties of which can 



be solved analytically. These ellipsoids were determined 
as best fits to the critical curves and caustics of the main 
galaxy haloes. The substructure population of each halo 
was taken directly from the simulations and superimposed 
on the analytic density fields of these isothermal ellipsoids. 
Multiple images of mock background quasars were recov- 



^ Defined as the radius enclosing a mean density of 200 times the 
critical density of the Universe. 



ered through an image-finding procedure combining the 
Newton-Raphson and Triangulation met hods (for a mo re 
detailed description of this procedure, see IXu et al.1 Hooi)). 
The cusp-caustic relations and their violation probabilities 
were calculated for the simulated lensing-galaxy haloes and 
compared to those derived from observations. In this paper 
we extend this method to study additional contributions to 
the lensing potential. 



2.3 Semi-analytic galaxies 

The Aquarius simulations provide us with particle distribu- 
tions of the main dark matter haloes and an associated pop- 
ulation of subhaloes. To investigate the effect of baryons in 
substructures, we pos tprocess Aquarius with t he GALFORM 
serni-analytic model (|Cole et al .11 1994 l2000l : I Bower et al.l 
Hoofi). We determine the dark matter me rger history of eac h 
subhalo directly from the simulations (|Hellv et al.l 1200^ ). 
From each of these 'merger trees', GALFORM computes the 
history of gas accretion, star formation and feedback in the 
corresponding galaxy. He re, as a fiduc i al mo del, we adopt 
the parameter values of iBower et al.l (|2006l ). which have 
been shown to be consistent wit h a number of obs erva- 
tional constraints. As discussed bv lCooper et al.l (|2010l '). the 
I Bower et al.l (|2006l ) parameters also agree well with the MW 
and M31 satellite galaxy LFs, provided that the value of the 
parameter Vcut is set to ~ 30kms~^ (a value consistent with 
recent estimates from numerical simulation fl) , in pr eference 
to the value of 50kms~^ used bv lBower et al.l(|2006l ). Adis- 



crepancy with observations is only apparent at 



> 



5, 



where constraints on the observed LF are weak. We stress 
that GALFORM is only used here to determine which dark 

^ Vcut is defined as the halo circular velocity below which a 
photoionizing background strongly suppresses further baryon ac- 
cretion and cooling. Rcionization is assum ed to occur instanta- 
neously at z = 6 in the lBower et al.l | |200(I) model. 
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matter substructures host stars, and to assign reasonable 
mass-to-light ratios to those substructures. For this purpose, 
a fiducial model producing a satellite LF consist ent with the 
availa ble data is sufficient. Our choice of the i Bower et al.l 
l|2006t ) p arameters with VcMt = 30kms~^ is identical to that 
used by ICooper et all ||2010D as a fiducial model in their 
study of galactic stellar haloes. Our findings concerning the 
lensing effect of substructures are not strongly sensitive to 
further changes in this model that could be made to match 
other properties of the satellite population. 

Populating the resolved subhaloes of Aquarius with our 
fiducial GALFORM model results in ~ 100 — 200 satellite 
galaxies within the virial radius of each Aquarius halo. As an 
example. Fig. [T] presents the properties of the semi-analytic 
galaxies (shown as blue symbols) in the halo Aq-A-2. The 
combined mass, Mtotai, of a satellite galaxy and its host 
subhalo ranges from ~ lO^h'^MQ to ~ 10^°/i"^Mq with 
a median of '-^ 3 x lO^/i~^M0. The baryonic mass Mbaryon 
spans more than six decades, which is consistent with the 
wide range of M W satellite luminos i ties (from ~ lOOL, 
to ~ lO^I/a, e .g. 



Strigari et al] l2008l : IKoposov etHllioq 



IStringer et al.ll2009l ). Their radial distribution has a median 
halo-centric distance of 70h~^ kpc. It is rare to find a satel- 
lite galaxy projected within the Einstein radius of the main 
halo (2 ~ 4/1"^ kpc). 



We require a structural model for each satellite galaxy 
in its host subhalo. Although GALFORM predicts the mor- 
phology and size of each simulated galaxy, we prefer to adopt 
a simpler and more conservative model that maximises the 
lensing effect. We describe all satellite galaxies using a sin- 
gular isothermal sphere (SIS) density profile. This profile, 
of the form p{r) oc r^^, leads to a constant defiection an- 
gle. The SIS profile of each galaxy terminates at a radius rt, 
which we set to be the half-mass radius rn of the galaxy's 
host dark matter subhalo (resulting in a median truncation 
radius of ~ 200 pc). The total stellar mass of the galaxy 
predicted by GALFORM is enclosed within rt. Under these 
assumptions, the median Einstein radius is ~ 0.006". Fur- 
ther increasing the concentration of our stellar components 
by setting rt = 0.3 rn (median rt ~ 60 pc) makes little 
difference to the final cusp-caustic violation probability (as 
discussed below, this is sensitive to the spatial distribution of 
the substructures in the main halo - see Fig. [TJ . All satellites 
have high mass-to-light ratios such that adiabatic contrac- 
tion due to the baryonic component can be neglected. The 
full hierarchy of sub-subhaloes within subhaloes is included 
in the lensing calculation, as in our earlier dark matter-only 
study. Galaxies hosted by these sub-subhaloes are also rep- 
resented by SIS profiles in this updated calculation. 

To determine the new cusp-caustic violations that re- 
sult from including baryons in dark substructures, the semi- 
analytic galaxies are added to the dark matter density fields 
from Aquarius. For the dark matter component of each sub- 
halo, the particle distribution from the simulation is used 
directly, as in our earlier study; we only assume SIS profiles 
for the stars. The new density fields are processed through 
our PM code and the image-finding routine to recreate lens- 
ing signals of background quasars. 



2.4 The cusp-caustic violation 

In any smooth lensing potential an asymptotic m agnification 
relation (the 'cusp-caustic re l ation'; [B landford fc Naravan 
19861 : ISchneider fc Weisslll992l : IZakharovlil995l : iKeeton et al 



20031) wiU hold: 



Reus 



MA + ^J■B + fj,c\ 



0, 



(1) 



\fJ.A\ + \f^B\ + ImcI 

where denotes the magnifications of the three closest im- 
ages A, B and C of a background source, which is located 
near a cusp of the central caustic. In this case, the total ab- 
solute magnification \fiA \ + Ia^bI -I- |mc| of the trip le images 
goes t o infinity. Following the same procedure as IXu et al.l 
|2009!), the lensing haloes are located at redshift 0.6 and 
the background quasars at redshift 3.0. 'Cusp sources' are 
defined as sources (quasars) that are near the cusps of the 
central caustic and have image opening angles, of their 
close triple images smaller than 90°. Any cusp source with 
^cusp 0.187 is defined as a cusp-violation 0.187 is 

the observed value for B1422-I-231, which shows the small- 
est violation (smallest -Rcusp value) among the five observed 
cusp lenses. In each of the three independent projections 
from each Aquarius halo, approximately 10000 - 20000 cusp 
sources are generated and used to calculate the -Rcusp distri- 
bution and the violation probability. 

Our earlier work jXu et al.i |200i) found an average 
cusp- violation rate of ~ 10% over all three independent pro- 
jections of each of the six haloe^, which leads to a proba- 
bility < 0.01 of observing 3 out of 5 cusp lenses violating 
the cusp-caustic relation due to substructures. This result is 
listed as case (a) in Table [T] which summarises the effects of 
additional contributions to the lensing potential. 

Case (b) in Table [1] adds the semi-analytic galaxies de- 
scribed above to the dark matter substructures. The change 
in the cusp-caustic violation is small. The reason for this 
is clear from Fig. [2] (left), which shows the critical curves 
superimposed on a contour map of the substructure sur- 
face mass fraction, in the halo Aq-E-2 (a single projection 
is shown as an example; our results are averaged over many 
projections). Th is figure can be compared with figure 11 in 
IXu et all (|2009l '): the addition of satellite galaxies alters the 
density profile of individual subhaloes, but cannot alter the 
abundance of substructures around the tangential critical 
curve, where the cusp-caustic relation is examined. There- 
fore, the baryonic 'boost' from the small number of galaxy- 
hosting subhaloes contributing to the overall cusp-caustic 
violations is marginal. 

As shown in Table 3 of IXu et al.l (|2009l '). the violation 
probability varies from halo to halo, and from projection to 
projection within the same halo. The lensing cross-section 
for cusp-caustic violation is significant only where massive 
substructures are projected onto the central region. The 
presence or absence of these strong perturbations to the lens- 
ing potential is highly stochastic, because massive subhaloes 

There were four projections in total with naked cusps (see 
IXu et aT]|2009l ) present in the central caustic due to the large core 
sizes of the imposed Hernquist profiles. The exclusion of these four 
cases in the final calculation resulted in a mean cusp-violation 
rate of 6.4%. Here in order to avoid these "naked cusps" , the core 
size is artificially set to be 0.05" for all six haloes, leading to the 
higher cusp-violation rate of ^ 10% as shown in Table[T]case (a). 
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Table 1. The average substructure surface mass fractions and the cusp-caustic violation probabihties for scenarios in which additional 
contributions to the lensing potential are considered. Col. 2: /sub.annu is the mean surface mass fraction of substructures within an 
0.1"-annulus around the tangential critical curve. Col. 3: P(-Rcusp ^ 0.187) is the mean cusp-violation rate averaged over all projections 
from six haloes, under which Col. 4, P315, gives the probability of observing 3 out of 5 cusp lenses violating the cusp-caustic relation due 
to the presence of substructures. 



Case 


/sub.annu 


P(Rcusp > 0.187) 




(a), dark matter subhaloes only 


0.20% 


10.1% 


0.8% 


(b), (a) -f semi-analytic galaxies 


0.21% 


10.2% 


0.9% 


(c), (b) -I- Milky- Way globular clusters 


0.23% 


11.0% 


1.1% 




Figure 2. Contour maps of the substructure surface mass fraction of the halo Aq-E-2, in X -projection. Left: semi-analytic galaxies are 
added to the dark matter subhalo population. The mean surface mass fraction, /gut, annu of the substructures within the 0.1"-annulus 
(indicated by the dashed lines) around the tangential critical curve is ^0.16%. Right: semi-analytic galaxies and Milky-Way globular 
clusters are added to the dark subhalo population, fsuh annu ~ 0.19%. Globular clusters are more centrally distributed in the projected 
central region than satellite galaxies. The blue circle indicates an example of the small-scale wiggles induced by globular clusters. 



are intrinsically rare and are typically found in the outer re- 
gions of their hosts. It is therefore possible in principle that 
a limited sample of haloes and projections could lead to an 
underestimate (or overestimate) of the violation probability. 
However, the mass function and radial distribution of sub- 
haloes is found to be very similar among all six Aquarius 
haloes (jSpringel et al.|[2008 ) . In these respects the Aquarius 
haloes are themselves repr esentative of the population from 
which they were selected (iBovl an-Kolchin et aLll2009l ). Our 
sample of 18 projections (combining three projections in six 
similar haloes) is therefore unlikely to significantly underes- 
timate the violation probability. 



3 MILKY- WAY GLOBULAR CLUSTERS AND 
THEIR CUSP-CAUSTIC VIOLATIONS 

Scenarios for the formation of globular clusters (GCs) are 
numerous and highly uncertain. GCs are known to have 
low mass-to-light ratios (i.e. to be baryon-dominated) at the 
present day (e.g. lMarastonll2005h . although a past associa- 



tion between individual GCs and dark matter haloes is not 
strongly ruled out. It may be that some GCs are the relics 
of the very fi rst generations of star formation in the early 
Universe (e.g. lCao et"alll2009l : [Criffen et al.ll2009h . It is of- 
ten speculated that a number of GCs (most notably uj Cen- 
tauri) should be identified with the stripped nuclei of dwarf 
satellite galaxies. This suggests that the definitions of these 
two classes of object might be ambiguous. At the present 
day GCs are strongly concentrated around central galaxies 
(more so than the overall subhalo populations in simulated 
CDM haloes) and their survival is known to be subject to 
many factors, including evaporation and tidal disruption. 
GCs are likely to be present in lensing galaxies in large 
numbers and may perturb the gravitational potentials in the 
inner regions, causing cusp-caustic violat ions. A simple es- 
timate of their contribution was made bv lMao &: Schneideij 
(1998), who concluded that a surface density fluctuation of 
a few per cent {5k ~ 0.01) from GCs would be enough to 
cause the observed flux-ratio anomaly in B1422-I-231. This 
conclusion has been re-examined in this work. 
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Here we adopt an empirical approach to the effects of 
GCs on the lens ing pot e ntial. We use the catalogue of Milky- 
Way GCs from iHarrid l|l996h . which provides their spatial 
distribution, V-band luminosities Lv and half-mass radii rh. 
Although the Milky- Way GC distribution is slightly flat- 
tened within the central ~ 10/i~^ kpc, the choice of projec- 
tion does not affect our results. 

It is interesting to ask whether a proportion of the 
'satellite galaxies' in Aquarius should in fact be identified 
with a population of 'primordial' galaxy-like objects or with 
the 'cores' of galaxies that have been heavily stripped. We 
have already included satellites in our calculation, so our 
approach to GCs risks double-counting some objects if ei- 
ther of these cases is true. However, a detailed investiga- 
tion of this issue is beyond the scope of this paper, and we 
will assume that most GCs are not already represented as 
'satellites' in our semi-analytic model. As we state above, 
the LF of bright satellite galaxies in our model matches the 
shape of the Milky- Way satellite LF. This observed LF does 
not include the many equally bright but structurally distinct 
Milky- Way GCs: this in turn suggests that these bright GCs 
are not represented by some of the existing 'satellites' in our 
model, stripped or otherwise. Fainter than Mv ~ —5 the 
distinction between galaxies and clusters is much less cer- 
tain and the interpretation of the current data is not at all 
clear. However, these low-mass objects are not significant 
for lensing. 

The masses of our empiric al Milky- Way G Cs are ob- 
tained assuming M/Lv = 3 fe.g. lMarastonll2005l ). We adopt 
a SIS profile for each GC to maximise our estimate of its 
lensing effect. We choose the SIS truncation radius rt = 2rh 
such that the model half-mass radius is equal to the GC 
half-light radius (Msis(< r) oc r). The median Einstein ra- 
dius is ~ 0.003" under these assumptions. Fig. [T] summarises 
our empirical GC population, showing the masses and radial 
distributions of the Milky- Way GCs (red symbols). Their 
dynamical masses range from hundreds to millions of solar 
masses (lower than that of the majority of satellite galaxies). 
The radial distribution of GCs is much more concentrated 
than satellite galaxies: the number of GCs projected at the 
Einstein radius of the main halo is about 1 ~ 2 times the 
number of the dark matter subhaloes at that radius. 

Comparing the right and left panels in Fig. [21 more sub- 
structures (GCs) in this case are projected in the vicinity of 
the tangential critical curves, inducing more smaller-scale 
wiggles; one example is indicated by the blue circle. Case 
(c) in Table [1] summarises our lensing results after includ- 
ing the empirical GC population. The mean cusp-violation 
rate averaged over all projections from six Aquarius haloes 
has only marginally increased from ~ 10% to ~ 11%, and 
the probability to observe 3 out of 5 cusp lenses violating 
the cusp-caustic relation due to substructures has increased 
from 0.9% to 1.1%. Although GCs are numerous and more 
concentrated than the subhalo population, the small mass 
and size of a typical GC limits their effect on the lensing 
potential. The predicted violation rate is still far below that 
observed, so the inclusion of GCs still does not solve the 
apparent lack of substructures suggested by observations of 
lensing flux anomalies. 



4 AQUARIUS STREAMS AND THEIR 
CUSP-CAUSTIC VIOLATIONS 

In recent years, streams and other stellar overdensities 
have been dete c ted in the halo of the Milk y Way (e.g. 
iBelokurov et all l2007l : iBelokurov et al.l 120071') and other 
nearby galaxies, most notably M31 fe.g. iMcConnachie et al.l 
I2OO9I ). Cosmological A'^-body simulations have demonstrated 
that a wealth of similar structures consistent with these 
observations are produced by the tidal st ripping of satel- 
lite galaxies in typical Ga lactic haloes (e.g. iHelmi fc White! 
I1999I : ICooper et al.llioTol '). As stars in satellites are likely 
to be deeply embedded in their own host dark matter 
(sub)haloes at the time of accretion, each stellar stream in 
these models is naturally associated with a more massive 
stream of dark matter. These 'streams' exhibit a variety of 
complex morphologies (e.g. Cooper et al 2009), very few 
of which remain 'coherent' in configuration space for long. 
However, their presence might induce irregularities in the 
surface density of the dark matter halo and cause lensing 
flux anomalies. 

We make an estimate of the magnitude of this effect 
with the level-2 Aquarius simulations. We define streams 
as dark matter particles that were once bound to substruc- 
tures crossing the virial radius of the main halo, but which 
at z — arc bound to the main halo itself (i.e. are no 
longer bound to any substructure). In the previous sections 
and in IXu et al.l (|2009|) , all these particles were treated as 
the smooth component of the main halo (included in the 
smooth ellipsoidal fit). The streams from different infalling 
substructures are phase-mixed to different degrees. Those 
most strongly mixed are smoothly distributed as ellipsoids 
concentric with the main halo, but a small quantity of satel- 
lite (dark matter) debris remains in the form of coherent 
features in configuration space (canonical 'streams'). These 
features introduce the most significant surface density irreg- 
ularities. 

Here we select the most massive spatially coherent 
streams that are projected within the central ~ Wh~^ kpc 
region of each halo. We do so by visually inspecting the den- 
sity fields of all massive streams (Mdm > 5 x W^h~^ Mq) for 
regions of high overdensity projected in the vicinity of the 
critical curve. We have calculated the lensing effect of each 
individual candidate stream by superimposing its full parti- 
cle distribution onto our smooth isothermal ellipsoid model 
for the main halo, as we did for the bound substructures 
described previously. In practice the overlap of many such 
streams at the centre of a lensing galajcy decreases the den- 
sity contrast of any single stream; by isolating these streams 
from the smoothly distributed dark matter, we are likely to 
overestimate their overall contribution. 

Fig. [3] shows one example of the massive coherent 
streams meeting our criteria. The surface density contours 
are overlaid by the lensing critical curves. The mean sur- 
face mass fraction, /stm.annu, of the stream within the 0.1"- 
annulus around the tangential critical curve, varies from 
0.01% to < 1% depending on projection. The cusp-violation 
rates in all cases are found to be extremely low, P(i?cusp ^ 
0.187) <C 1%, even for projections with /stm.annu as high 
as those from compact substructures. The typical width of 
these 'dens e' spatially coherent s treams is of the order of 
~ kpc fe.g. iHelmi fc White|[l999l ): this is about the size of 
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Figure 3. An example of one spatially coherent stream in the halo Aq-A-2. The top two and bottom left panels show the stream on 
large scales in three independent projections. The bottom right panel presents the contour map of the surface mass density fraction of 
the stream on smaller scales in Y-projection, with critical curves overlaid; /stm.annu ~ 0.45%. 



the strong lensing domain and hence these streams are inca- 
pable of inducing small-scale wiggles in the tangential criti- 
cal curves. On the other hand, thin streams from lower mass 
progenitors are not sufficiently dense. We conclude that the 
presence of dark matter streams does not significantly boost 
the cusp-caustic violations. 



5 DISCUSSION AND CONCLUSIONS 

We briefly discuss a few issues which may affect our conclu- 
sions. Firstly, the observed sample is small, consisting of only 
five cusp lenses, among which are the thre e systems detected 
in the radio band b y the CLASS survey (|Mvers et al.ll2003l : 
iBrowne et alil2003l ). This survey shows an ano malously high 
incidence of secondary len ses (bright satellites) (| Bryan et al.l 
l2008l : I Jackson et al.|[201ol ). 

A larger sample of galaxy-scale lenses from fu- 
ture instruments will provide better statistical constraints 
on substructure abundance in galaxy-sized CDM haloes 
(jVeeetti fc KoopmansI I2OO9I I. New lenses are likely to be 
discovered in optical surveys (such as PanSTARRS and 
LSST), although radio and mid-infrared follow-up will be 



necessary to distinguish anomalies attributable to substruc- 
ture from microlensing events (e.g. IWambsganss et al. II1990I : 
IWitt et al. I I 19951 : iPoolev et al. l l2007l : lEigenbrod et al.ll20oi ). 
High-resolution radio surveys (to be carried out by LOFAR 

and SKA, for example) will be a productive way to discover 
more of these systems. 

Secondly, most lensing galaxies are massive ellip- 
tical galaxies due to their large lensing cross-sections 
iTurner et al.lll984 ). Compared with their spiral and irregu- 
lar count erparts, massive ellipticals are known to host m ore 
GCs (e. g. iForte et"aDll982l : iHarri J 199ll . lT99^ : l WestI 199^ 1 . In 
this work, we consider the lensing effect of GCs by imposing 
the observed Milky- Way GC population. If the lens galax- 
ies are much more massive than the Aquarius haloes, this 
assumption may underestimate the appropriate GC abun- 
dance. However, even with an increased number of GCs typi- 
cal of massive ellipticals, we would not expect major changes 
to our conclusions. 

Thirdly, the semi-analytic galaxies that we have in- 
cluded in the lensing simulation are all hosted by dark mat- 
ter subhaloes above the resolution limit of the Aquarius sim- 
ulations. However, semi-analytic models also follow the at- 
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tached stellar components even after their host dark mat- 
ter (sub)haloes are tidally stripped below 20 particles 
IQ^H'^Mq). The baryonic masses of these 'sub-resolution' 
stellar components are in the same range as those of satellite 
galaxies, and each is associated with a bound dark struc- 
ture no more massive than lQ^h~^ Mq. In GALFORM their 
survival depends on an estimate of their merger timescales, 
which is subject to the treatment of dynamical friction. The 
radial distribution of the centres (most-bound particles) of 
these objects is more concentrated than that of their coun- 
terparts with resolved dark haloes. Nevertheless, including 
these sub-resolution 'orphan galaxies' only marginally in- 
creases the cusp- violation probability under the assumption 
of singular isothermal spheres for their density profiles. 

Finally, when including semi-analytic galaxies to 
account for the lensing effect of baryonic substructures, we 
have assumed that the overall galactic potential and the 
dynamics of subhaloes would not be significantly modified 
if baryons were added to our dark matter-only simulations. 
However, a concentration of baryons at the centre of a 
subhalo may increase its resilience to tidal stripping. The 
abundance of subhaloes in the inner region of the main 
halo may increase as a result of this effect. Our model 
does not account for this 'baryon-enhanced' survival of 
a small number of massive objects, although we do not 
consider this omission to be significant. There are likely 
to be few such objects, as a subhalo must typically lose 
more than ~ 90% of its mass before the structure of its 
central (stellar) c ore influences further stripping (see e.g. 
iPefiarrubia et al. lioos). This naturally limits the number 
of 'marginally destroyed' cases that would change in a fully 
self-consistent model. Another effect of including baryons 
would be to increase the concentration of the main halo in 
our simulation, resulting in a stronger tidal field. Ifowever, 
studies of disc-subhalo interactions using the Aquarius 
haloes (Lowing et al. in prep) support our assumption 
that the subhalo population (in the regime relevant to 
lensing) is only marginally affected by the presence of a 
realis t ic stellar disc. D isc shocking (e.g. iKazantzidis et aU 
I2OO9I : I D'Onghia et al.l 2009) may also reduce the inner 
substructure abundance, which would further aggravate the 
problem of explaining the observed lensing anomalies. 

et al.l ([200^) , the dark matter-only Aquarius sim- 
ulations were used to study the relation between dark matter 
substructure abundance and lensing flux-ratio anomalies, in 
particular the cusp-caustic violations observed for multiply 
imaged quasars. We found that the dark substructures in- 
trinsic to a typical galaxy-scale lensing halo were not suffi- 
cient to explain the observed frequency of cusp-caustic vio- 
lations. In this work, we have considered whether this expec- 
tation changes when satellite galaxies in subhaloes are taken 
into account, or when the lensing halo is assumed to contain 
a MW-like globular cluster population that is not associated 
with surviving dark matter substructures. We have also con- 
sidered streams of dark matter identified in the Aquarius 
simulations, in order to estimate the lensing effect of irregu- 
larities in the halo itself. We conclude that the abundance of 
intrinsic substructures, dark or bright, bound or diffuse, can- 
not fully account for the observed cusp-violation frequency. 
Taken at face value, this lack of substructure suggests a se- 
rious problem for the CDM model. Warm dark matter mod- 



els, which could reduce the satellite abundance and may 
help to bring the dwarf galaxy LF into agreement with ob- 
servations without invoking strong feedback or photoioniza- 
tion effects, (e.g. ,Sawala et al . 2010) would only make this 
problem worse. However, it is possible that the observed 
frequency of flux anomalies are strongly biased by the small 
number statistics of CLASS. 

Previous studies have shown that intergalactic haloes 
{M < IQ^'^Mq) projected along the line of sight can cause 
surface density fluctuations at the level of 1-10 per cent, 
making them a probabl e source of lensing flux anomalies 
(Chen et al. 2003; Wam bsganss et al.ll2005l : lMetca lf 2005a Q; 
Miranda & Maccio 2007). In future work, we will examine 
the contribution of this large-scale structure to observations 
of the cusp-caustic violation rate, using high-resolution cos- 
mological volume simulations, which self-consistently model 
the halo mass function and clu stering along the line of sight 
(e.g. iPuchwein fc Hilbert|[2009l ). 
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